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ABSTRACT 

The currently most viable methods to estimate magnetic field strengths in molecular cloud cores 
are Zeeman measurements and the Chandrasekhar-Fermi (CF) method. The CF-method estimates 
magnetic field strengths from polarimetry and relies on equipartition between the turbulent kinetic and 
turbulent magnetic energy in the observed region. Thus, its application to objects not dominated by 
hydromagnetic turbulence is questionable. We calibrate the CF-method and Zeeman measurements 
against numerical models of self-gravitating molecular cloud cores. We find - in agreement with 
previous results - that the CF-method is on the average accurate up to a factor of 2, while a single 
estimate can be off by a factor of 7. The CF-method is surprisingly robust even when applied to regions 
not dominated by hydromagnetic turbulence and in the face of small-number statistics. Zeeman 
measurements can systematically overestimate the field strength in high-density regions, if the field 
is predominantly oriented along the line of sight, while they tend to underestimate the field in other 
cases. 



Subject headings: polarization — MHD 
— ISM:globules 

1. MOTIVATION 

Magnetic fields are omnipresent in the interstellar 
medium (ISM). Their importance depends on their 
strength. The large-scale Galactic magnetic field con- 
tributes to the scale height of the Galactic disk, and 
it confines cosmic rays to the Galaxy, thus influenc- 
ing the chemistry in molecular clouds. Molecular cloud 
cores - thought to be gravitationally unstable with- 
out an additional support mechanism - have been sup- 
posed to be stabilized by magnetic fields, leading to the 
mod el of magnetical ly mediated (low mass) star forma- 
tion llShu et al.ll987|) . An alternative - althoug h not nec- 
essarily exclusively so IjLin fc Nakamurall2004|) - model 
pictures magnetic fields as dynamically unimportant, 
and star formation being controlled mainly by turbulence 
on dynamical time scales (e.g. Larson 1981; Elmegreen 
2000; Hartmann et al. 2001; Mac Low k. Klessen 2004). 

Magnetic field strength estimates in molecular cloud 
cores thus might provide a way to distinguish between 
those two models. Such measurements rely mostly on 
the Zeeman effect (e.g. Crutcher 1999, Bourke et al. 
2001, Sarma et al. 2002) resulting in an estimate of 
the density-weighted, direction-sensitive integrated mag- 
netic field strength along the line of sight, and on the 
Chandrasekhar-Fermi method ( Chan drasckhar fc Fcrnij 
1195,'^ . The CF-method relates the line of sight velocity 
dispersion to the plane of sky polarization angle disper- 
sion under the assumption of an isotropically turbulent 
medium whose turbulent kinetic and turbulent magnetic 
energy components are in equipartition, yielding an esti- 
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In its original form as given by iChandrasekhar fc Fermi 
l|1953l) . equation (Q) assumes small angles around the (re- 
quired) mean field direction. Calibrations of the method 
with the help of numerical simulations of magnetized tur- 
bulent molecular cloud regions show that the CF-method 
on average overestimates the "true" field strength by a 
factor of approximately 2 (Ostriker et al. 2001; Padoan 
et al. 2001; Heitsch et al. 2001b, in the following 
HZMLN). Limited observational resolution (mimicked 
by smoothing the simulated polarization maps with a 
Gaussian) leads to a further systematic overestimate 
(HZMLN). While the meth od is reasonably robu st when 
selecting for small angles l|Ostriker et al.ll2001|) . it fails 
utterly for vanishing mean fields in the plane of sky, as 
expected (HZMLN). The latter authors suggested a mod- 
ification of the CF-method which estimates not the mean 
field (B), but the rms field {B^^^- 

The CF-method complements Zeeman estimates to 
some extent, since it infers the mean field strength in the 
plane of sky. The hope is to actually retrieve information 
about the 3D structure of the magnetic field (Lai et al] 
I2003J) . Many of the more recent applications of the CF- 
method refer to the dense phases of the ISM, mostly to 
molecular cloud cores or Bok globules (e.g. Lai et al. 
2001; Henning et al. 2001; Lai et al. 2002; Matthews 
et al. 2002; Lai et al. 200 3, Wolf et al. 2003; Crutcher 
et al. 2004). iVallee et all (|2003) decided not to use the 
CF-method because of the uncertainty introduced by the 
line-of-sight averaging of the polarization pattern. Since 

* There is a difference between the mean field with respect to di- 
rection and the mean field with respect to orientation. The former 
refers to a local field perturbation (which could be perpendicular 
to the mean field w.r.t. orientation), while the latter is an averaged 
mean field. It is the latter we are referring to. 
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the method rehes on equipartition between the turbu- 
lent kinetic and magnetic energy, its application to other 
regimes than those dominated by hydromagnetic turbu- 
lence might raise some questions that we will discuss in 

The goal of this study is twofold: First we aim at es- 
tablishing how reliable the CF-estimates are for regions 
which are not dominated by hydromagnetic turbulence. 
We calibrate the method with the help of cores in a 
high resohition model of self-gravitating magnetized tur- 
bulence (Li ct al. 2004). We find that the CF- method 
returns on average reliable estimates (within a factor of 
2) as long as the full angle information is used. Single 
estimates can be off by a factor of 7. The method is sur- 
prisingly robust even when it is applied to rotating and 
self-gravitating cores. Despite its robustness, the method 
is subject to various systematic effects. 

Second, we use the model cores to infer Zeeman mea- 
surements of the field strengths and to compare those to 
the CF-estimates. While Zeeman measurements are ex- 
pected to underestimate the magnetic field in a turbulent 
region because of cancellation along the line of sight, we 
find that they can systematically overestimate the field in 
high-density regions with ordered fields oriented parallel 
to the line of sight. 

2. THE CF-METHOD IN CORES: SYSTEMATIC EFFECTS 

The CF-method in its original form rests on four 
assumptions: (a) The field perturbations are caused 
by hydromagnetic turbulence, (b) The turbulent ki- 
netic and turbulent magnetic energy are in equiparti- 
tion (This is guarantee d for A lfvenic perturbations, see 
e.g. iZweibel k McKed IjlQQSft .V (c) The turbulence is 
isotropic, (d) The field is perfectly frozen to the (to- 
tal) gas density. If one or more of these assumptions are 
violated, systematic effects will be introduced in the CF- 
estimates (see e.g. Zweibel (1990), Myers & Coodmann 
(1991), Zweibel (1996) and Crutcher et al. (2004) for 
a discussion). In the context of protostellar cores, two 
effects would be expected to play a major role: 

(i) The field perturbations are caused by other agents 
than turbulence, namely self-gravitation and/or rotation. 
First, we could imagine a core within the picture of mag- 
netically dominated star formation: The core is perme- 
ated by a strong field, rendering the system subcritical 
(e.g. [Mouschovias & Soitzcr (1976)). Turbulence is as- 
sumed to be negligible. The field would be expected 
to evolve to an hourglass-like structure, which the CF- 
method would in turn interpret as an effect of turbulence, 
thus underestimating the field strength. 

Second, let us consider a fiattened rotating core seen 
face on such, that the rotation axis is more or less parallel 
to the line of sight. The mean magnetic field vector is 
oriented parallel to the line of sight (This is the situation 
envisaged in the concept of magnetic brakingof cores, 
see e.g. [Liist & Schhitcr (1955); Mcstcl (195^)). The 
field lines would wind up, i.e. in projection they would 
describe spirals. These would result in large polarization 
angle dispersions not necessarily accompanied by a large 
velocity dispersion along the line of sight. 

(ii) For anisotro pic turbulence (see e.g. 
iCho fc LazarianI ()200.lD 'l. relating the line of sight 
velocity dispersion to the polarization angle dispersion 
can introduce a systematic effect in the CF-estimate. 



However, if the magnetic field is predominantly in the 
plane of sky, the CF-method could still yield reasonable 
results despite severe anisotropy. 

3. MODELS AND METHODS 

The m odel this study rests on is described in detail by 
iLi et alJ (|2(304). It mimics a fraction of an isothermal 
molecular cloud region (periodic boundary conditions), 
which is dominated by supersonic turbulence driven on 
the largest scales at Mach 10 such that the driving en- 
ergy inpu t is kep t constant (for a detailed description see 
iMac Low! ()1999j) V Once an equilibrium state between 
driving energy input and dissipation has been reached, 
self-gravity is switched on. The shock-generated fila- 
ments start to fragment and form cores, the objects 
of this study. The initial ratio of thermal to magnetic 
pressure is /3 = 0.9. For a detailed discussion of this 
type of simulation an d the choice of parameters see e.g. 
iHeitsch et all (|2001aD . 

Cores selected for this study have to fulfill several cri- 
teria: (a) They are gravitationally bound, (b) They 
are resolved accordin g to the so-called Jeans-criterion 
(|Trueloveet al.llT997l) . given by 



A., 




> AAx, 



(2) 



where Ax is the grid spacing and Cs is the sound speed, 
(c) Regions of surface densities of 10% of the peak sur- 
face density or above are used for analysis. Although 
this (excitation-) threshold is chosen somewhat arbitrar- 
ily (albeit guided by observational maps), other choices 
(within reasonable limits) did not affect the results, (d) 
Selected cores consist of more than 1000 cells in order to 
sample their internal structure sufficiently. The sample 
cores are listed in Table Q] The fourth column in the 
table denotes unresolved cores according to the Jeans 
criterion. These cores are not used for the subsequent 
analysis except to demonstrate resolution effects. Note 
that p/pj gives a resolution criterion, while M/Mj deter- 
mines whether a core is gravitationally unstable. Here, 
Mj is the thermal Jeans mass. 

Summarizing Tabled the model cores are marginally 
gravitationally bound, their internal velocities are (in 
most cases) subsonic {M « 0.5), and the magnetic pres- 
sure is on the order of the thermal pressure (with a sub- 
stantial scatter). 

Once cores have been identified over the run of the 
simulation, we construct surface density maps and po- 
larization maps following iZweibeJI ()1996|) via 

fiv) D2 , J cos" idy (3) 



P = Q + iU = N 



Bl + Bl 

where B^ and B^ give the field vectors in the plane of 
sky perpendicular to the line of sight and are taken di- 
rectly from the simulations. We integrate along the line 
of sight in y-direction. The function /(y) is a weighting 
function which accounts for the density, emissivity, and 
polarizing properties of the dust grains. As the simplest 
assumption, we take f{y) to be the local gas density nor- 
malized by the column density 1/A^ = / f{y) dy, where 
we take the opportunity to correct a typo in HZMLN. 
The factor 



cos^ 7 
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Table 1. Properties of sample cores 



core 


pIpj 


# cells 




M/Mj 


M 


/3 


Ma 


*/^// 


d04c04 


0.32 


5956 




1.4 


0.2 


2.2 


0.2 


0.50 


d05c02 


0.33 


5283 




1.3 


0.3 


1.9 


0.3 


0.67 


d06cll 


0.28 


3589 




1.6 


0.5 


3.7 


0.7 


0.83 


d06cl5 


0.22 


2928 




0.8 


0.5 


1.1 


0.4 


0.83 


d07c04 


5.68 


2167 


X 


5.4 


1.0 


2.0 


1.0 


1.00 


d07c05 


2.44 


1194 


X 


2.7 


0.7 


2.3 


0.8 


1.00 


d07cll 


0.55 


1574 




1.8 


0.6 


0.8 


0.4 


1.00 


dOScOl 


70.6 


2043 


X 


23 


1.4 


0.7 


0.8 


1.17 


d08c02 


57.6 


2056 


X 


21 


1.2 


1.9 


1.2 


1.17 


d08c09 


9.58 


1133 


X 


6.5 


1.0 


1.6 


0.9 


1.17 


d08cl7 


0.26 


2344 




2.2 


0.3 


6.8 


0.6 


1.17 


d08cl9 


0.26 


1354 




0.5 


0.3 


8.1 


0.6 


1.17 


d09cl4 


0.91 


1440 




1.8 


0.4 


1.5 


0.4 


1.33 


d09cl9 


0.35 


1130 




1.2 


0.2 


3.3 


0.3 


1.33 



Note. — First column: core label: dicj, where i denotes the 
time dump number and j specifies the core number at time dump 
i. Second column: maximum over Jeans density, derived from 
the stability criterion following ITruelove et al. 11997). Unresolved 
cores are marked with an X in the fourth column. Third column: 
number of cells within core. Fifth column: mass contained in core 
over Jeans mass, assuming a spherical gas distribution. Sixth col- 
umn: Mach number derived from (density-weighted) line-of-sight 
velocity dispersion. Seventh column: plasma beta /3 = 2pc^/_B^. 
Eighth column: Alfven Mach number Mj\ = i^rms/cA where the 
Alfven speed c^ = B/y/ATvp. Ninth column: dump time in units 
of the global free fall time Tf f = y^3n/(32Gp) ~ lO^yrs. 



accounts for suppression of polarization by the mag- 
netic f ield component along the l i ne-of- sight (e.g. 
iFiege fc Pudritz (2000), Padoan et aP ^Wlh)- The po- 
larized intensity is \P\ ~ \/Q^ + V^ , and the polarization 
angle is 

Q ~ — arctan -- . 5 

2 Q ^ ' 

Equation Q is an approximate solution of the full radia- 
tive transfer equation for the Stokes parameters, (Mar- 
tin 1974; Lee & Draine 1985), valid for small polariza- 
tion and low optical depth. At far infrared wavelengths, 
and at typical column densities for molecular clouds, 
the medium can safely be assumed to be optically thin 
(^Iildcbrand et al. 2000). We deliberately leave out any 
statement about the polarization degree as we are only 
interested in the polarization angles, but not in the ra- 
tio between polarized and continuum intensity. By this 
we intend to keep the argument as simple as possible, 
not embarking on a discussion o n the effic iency of grain 
alignment mechanisms (see e.g. iLazarianl II2003J) '). This 
simplification certainly limits our study to the investiga- 
tion of polarization angles only. To mimic limited tele- 
scope resolution, we smooth the polarization maps with 
a Gaussian filter as in HZMLN. 

The CF-method is applied to the "polarization maps" 
in the extended form given by 



(B)2 = 47r/9 



o{vf 



cr(tan6l)2 



(6) 



where we replace the small-angle approximation by the 
tan of the angle. The density is determined from the 
column densities of the polarization maps, assuming a 
spherical mass distribution, where the radius is given by 
the geometric mean o f the long and short axis of the core 
IjCrutcher et al.ll2004|) . 




Fig. 1. — 3D core projection. The core outline is given by iso- 
density surfaces. The straight line denotes the mean flux direction. 
Three representative field lines are outlined. 



To infer a "Zeeman" estimate from our cores, we inte- 
grate 



{Bzm)^N fiy)Bydy 



(7) 



along the line of sight. N and /(y) are the same as in 
equation (|3J|. 

4. RESULTS 

In order to establish the reliability of the CF-method 
for cores, we generated polarization maps after rotating 
the cores such that the line of sight is either parallel to 
the mean flux direction (the "worst case" for the CF- 
method, since then the necessary mean field component 
would be expected to vanish) or perpendicular to the 
mean flux direction ("best" case). An example of the 
three-dimensional core structure and the resulting polar- 
ization maps is given in Figures ^ and |21 The depicted 
core is a typical example in the sense that it combines 
very ordered with turbulent field structure. Figure |21 
displays the three polarization maps - again seen along 
the mean flux direction and perpend i cular to it - for the 
core shown in Figure 22 of iLi et all l)2004l) . The oblate 
core is permeated by field lines oriented in parallel to its 
rotation axis. The field lines are slowly wound up, result- 
ing in a spiral pattern when viewing the core along its 
rotation axis (left panel of Fig. O). For the lines of sight 
perpendicular to the rotation axis, the rotation pattern 
does not emerge (center and right panel). 

4.1. Performance of the CF-method in cores 

Figure 01 summarizes the reliability estimates for the 
CF-method applied to the resolved model cores (i.e. the 
ones without X's in Table^. For this plot, all three lines 
of sight were used. First, we note that using the small 
angle approximation (open symbols) leads to a system- 
atic overestimate of the field strength by a factor of more 
than 4. This is not surprising since we discard physically 
relevant information about the field structure (and thus 
strength) when selecting for small angles. However, we 
cannot restrict the line of sight velocity dispersion in the 
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Fig. 2. — (a) Polarization map looking along the mean flux direction. (b,c) Polarization map looking perpendicular to the mean flux 
direction denoted in Fig. One panel measures approximately 0.25pc across. The color bars give the column density in particles cm~'^. 
Every second polarization vector is shown (and used for the field estimate), resulting in a spatial "resolution" of 0.012pc. 
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Fig. 3. — (a) Polarization map looking along the mean flux direction for the core shown in Fig. 22 of lLi et al.l <2004l) . Note the spiral 
polarization pattern induced by the core's rotation. (b,c) Polarization map looking perpendicular to the mean flux direction. One panel 
measures approximately 0.25pc across. The color bars give the column density in particles cm~-^. Every second polarization vector is 
shown (and used for the field estimate), resulting in a spatial "resolution" of 0.012pc. 



same way. On the other hand, using equation JB)) (filled 
symbols) leads to a moderate overestimate of the field 
by a factor of approximate ly two^. This i s not incon- 
sistent with the findings of lOstriker et al.1 l)200H) . since 



they selected for their whole domain and thus their re- 
sults were not affected by small-number statistics. For 
large-number statistics, weak fields will lead to a wide an- 
gle dispersion, meaning that even when selecting, most of 
the angles are found at the largest possible values, thus 
indicating a weak field. 

The small number of strong-field cores is conspicuous. 
This is simply a consequence of low initial field strength, 
which renders the turbulent box supercritical by a fac- 
tor of 8. On the one hand, this is inconsistent with the 
observations, on the other hand it allows us to test the 
CF-method under severer conditions than it would prob- 
ably meet in observational data. 

The modified method to estimate the rms field in the 
core gives - within the scatter - the same estimated field 
(Fig. ISJ. For smooth fields the method reverts to the 



5 The selection for Ad < 85° 
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Fig. 4. — Reliability a = Bcf/{B) versus the mean field in 
the core {B) for all resolved cores and lines of sight. Open sym- 
bols: angular dispersion selected for A0 < 25°. Filled symbols: 
angular dispersion selected for A6 < 85° . The dashed lines denote 
the averages, whose values are given in the plot. No smoothing 
applied. Clearly, selecting for 25° leads to a systematic overesti- 
mate, while using the full information leads to an overestimate by 
approximately 2 (lower dashed line). 
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Fig. 5. — Reliability a = Brms/{B'^)^'^ versus the rms mean 
field in the core (i?-^)^" for all resolved cores and lines of sight. 
Open symbols: angular dispersion selected for A9 < 25°. Filled 
symbols: angular dispersion selected for A9 < 85°. The dashed 
lines denote the averages, whose values are given in the plot. No 
smoothing applied. 



original version, while for unordered fields, the method 
emphasizes the angular variations (see HZMLN for a dis- 
cussion) . The scatter of the estimates is slightly reduced. 
Note that on the whole the rms field is somewhat larger 
than the mean field, indicating a non-unifornr field com- 
ponent in the cores. 

Figures 0] and all subsequent ones show cores ex- 
tracted at different times during the simulation. Al- 
though we thus do not expect to be affected by selec- 
tion effects that could result from chosing a special time 
instant during the simulation, there is of course a se- 
lection effect in the sense that we extract cores at com- 
parable evolutionary stages, i.e. cores that are already 
self-gravitating, but still numerically resolved. As long 
as the cores are resolved, a correlation between reliabil- 
ity and global evolution time (last column in Table ^1 
cannot be seen. 

Applying equation @ and preventing the angular dis- 
persion from blowing up by selecting for angle differences 
A9 < 85° improves the reliability of the estimates in 
these models. 

4.2. Effects of limited observational and numerical 
resolution 

As HZMLN showed for extended regions, limited tele- 
scope resolution (finite beam width) leads to a systematic 
overestimate, because the field variations are averaged 
out, making the field look smoother, i.e. in terms of the 
method, stronger. We can observe a similar effect in Fig- 
ureini compared to Figure^ Degrading the resolution by 
roughly a factor of 3 leads to an increased overestimate. 
The estimates selected for A9 < 25° are less affected by 
the degraded resolution, and both selected sets of esti- 
mates agree for strong fields, as is to be expected. 

It might be interesting to note that the ambipolar dif- 
fusion scale 



L^^^h^^ 1.97 X 10^1 



S^ 



(8) 
is on the order of O.Olpc for typical core parameters. In 



equation IJSJ the field strength B, the neutral particle 
density ri„ and the velocity dispersion Vrms are given in 
cgs- units. The ionization fraction is denoted by Xi. 
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Fig. 6. — Reliability a = Bqp/{B) versus the mean field in 
the core (B) for all resolved cores and lines of sight at reduced 
resolution (map smoothed with a Gaussian of cr = 3pix). Open 
symbols: angular dispersion selected for AS < 25° . Filled symbols: 
angular dispersion selected for A6 < 85° . The dashed lines denote 
the averages, whose values are given in the plot. 

While Figure El demonstrated the effect of degraded 
telescope resolution, we exemplify the infiuence of under- 
resolved cores on the overall result in Figured The ef- 
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<B> [/xG] 

Fig. 7. — Reliability a = Bcf/{B) versus the mean field in the 
core (B) for all sample cores (including unresolved ones, see Tab. IT] 
and lines of sight. Open symbols: angular dispersion selected for 
A8 < 25°. Filled symbols: angular dispersion selected for A8 < 
85° . The dashed lines denote the averages, whose values are given 
in the plot. 

feet is obvious: The peak densities in those cores start 
to grow unhampered once they overstep the resolution 
criterion. The peak is concentrated in the central cell, so 
that the field is not a ble to follow the density evolution 
ll Heitsch et al.ll2(Tniaj) . Because of the over-large den- 
sities, the CF-method systematically overestimates the 
true field strength, a purely numerical effect in this case. 
Note that the overestimates correlate with the strongest 
fields and unresolved cores. 

The insufficient numerical resolution leads to a diffu- 
sion of the field out of the dense region, qualitatively 
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Fig. 8. — Reliability for the Zeeman estimates Bzm /{B) against 
the mean field in the core. Averages broken up according to line 
of sights are given in the plot, and los=l means integration along 
the mean flux direction (as in Fig. l2l. 



similar to ambipolar diffusion. While there is no ba- 
sis on which to draw any conclusions from these under- 
resolved cores (especially since the numerical diffusion 
cannot be controlled, and can behave quite differently 
from its physical counterpart), a similar effect would be 
expected in the observations if the density tracer is a 
neutral species. 

4.3. CF Method vs. Zeeman Effect 

If the magnetic field is not uniform, Zeeman measure- 
ments - realized in our models via equation lO - are ex- 
pected to systematically underestimate the field strength 
because of cancellation by field reversals along the line 
of sight. At first glance, this expectation is not met 
(Fig. inj. The Zeeman-effect seems to overestimate the 
field roughly by the same amount as the CF-method 
does. However, distinguishing between the lines of sight 
clarifies the picture. Zeeman estimates taken along the 
mean flux direction - denoted by "los=l" in Figure |H1 
- yield larger values than those taken perpendicularly to 
the mean flux direction ("los=2" and "los=3" ). Since the 
line of sight integration is density-weighted, fields at lo- 
cations with higher densities tend to be emphasized. Be- 
cause of flux- freezing, this generally means that there is a 
tendency to higher field strength estimates in cores with 
large density contrasts, leading to a systematic overesti- 
mate for such cores. Thus the estimates along the mean 
fiux direction are dominated by the strongly weighted 
central field values, while the estimates perpendicular to 
the mean flux direction "see" the perturbed component 
around it, leading to cancellation. This of course holds 
only as long as the emissivity is (approximately) propor- 
tional to the density, as assumed in equation (0). 

Comparing the CF and Zeeman estimates in Figure |5| 
the picture is somewhat less clear for the CF-method. 
There is a slight tendency to a larger scatter for es- 
timates along the mean flux direction (upper panel of 
Fig.O diamonds), as one would expect, since the method 
is expected to work most reliably for a strong mean field 
component in the plane of sky. Unfortunately, this is a 
reliability criterion gathered in "hindsight" , i.e. it is not 
possible to apply it to observational data. 



CQ 

V 



CQ 

V 



1 



1 


: <> 


■ 




♦ 


■ 




-i 


r 


o 


o«°° n 




O 

D 


• 


i 




: 0_ 




s^ ■ ^. 


-♦ 




Zccrnan 


^ 


■ ■ 

° 8 • 




■ 

O 


^■•■■■■ 


1 






° oo 




D 


« los=l 
■ los = 2 
• los = 3 


: 







20 



40 



60 



80 



100 



Fig. 9. — Reliability for the CF-estimates as in Figure |4] (upper 
panel) and for the "Zeeman" estimates Bzm/{B) (lower panel) 
against the mean field in the core. The dashed line gives the mean 
reliability. los=l means integration along the mean flux direction 
(as in Fig. |2J. Filled symbols denote estimates for which Bqp 
and Bzm agree within a factor of 3, open symbols stand for the 
remaining estimates. 



Instead, we might be tempted to argue that for cases 
where the CF-method and the Zeeman measurements 
yield approximately the same results, the field value 
should be well determined. Filled symbols in Figure El 
denote estimates for which Bcf and Bzm agree within 
a factor of three. The scatter is reduced for both the 
CF-method and the Zeeman- measurement, although we 
should note that the estimates could still be off by a 
factor of approximately five. It might be interesting to 
note that the overall scatter of the Zeeman measurements 
around the "true" field value is actually larger than that 
of the CF-estimates (Fig. |2J. 

4.4. Effects on derived criticality of cores 

Figure El demonstrates the scatter of the magnetic 
field estimat es in terms of the cri ticality parameter A 
as defined bv lCrutcher et al.l lj2004|) . eq.(l), namely 



A = 7.6 X 10 



-21 ^(H2) 

B '' 



(9) 



where the magnetic field B is given in /iG and the H2 
column density in cm^^. 

The tendency of the CF-method (squares in Fig.ll()|l to 
generally overestimate the field strength renders the core 
ensemble slightly more subcritical than expected from 
the numerical values (diamonds in Fig. ll()|l . There are a 
few subcritical numerical data points (diamonds). These 
are strongly elongated filaments with fields aligned in 
parallel, indicating field line stretching. Although the 
numerical values already show a considerable scatter be- 
cause of dynamical effects, the uncertainties in the CF- 
estimates as well as in the Zeeman-estimates emphasize 
that field estimates of single cores allow only a limited 
conclusion about the physical state of the core. 

5. DISCUSSION 

Zeeman measurements and the CF-method are cur- 
rently the most viable ways to estimate the magnetic 
field strength in the dense molecular medium. Since the 
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Fig. 10. — Column density A'^ against magnetic field strength 
B to estimate criticality of cores. The line gives A = 1 (see 
eq. [1] of Crutchcr ct al. (2004)). The symbols denote how N{B) 
is determined. Diamonds: directly from simulations, squares: CF- 
Method, circles: Zeeman effect. Only resolved cores are shown (see 
Tab.0, selected for A6 < 85° in the case of the CF-estimates. Dif- 
ferent symbols on constant lines of A^ belong to one line of sight 
of one core. The star denotes N{B) for the full simulation domain 
(i.e. 27 pc^). 



two methods are measuring different components of the 
field, they complement each other to some extent. The 
CF-method rests on the assumption of equipartition be- 
tween turbulent kinetic and turbulent magnetic energy 
in the medium, which is not necessarily guaranteed in 
regions dominated by e.g gravity or rotation, such as 
protostellar cores. 

We tested with the help of model cores formed in a nu- 
merical simulation of self-gravitating magnetized turbu- 
lence how accurately the CF-method and Zeeman mea- 
surements can estimate the magnetic field strength in 
such cores. We find 

1. that the CF-method is on the average surprisingly 
robust even in regions not dominated by turbu- 
lence. The previous overestimation factor of 2 is 
still valid approximately for the model cores on the 
average, ?/the full angle dispersion information is 
used. Single measurements can be off by a factor 
of up to 7 in both directions. 

2. that applying the small-angle approximation in the 
CF-method leads to a systematic overestimate of 
{B). 



3. that the "Zeeman" estimates of our cores can 
lead to an overestimate for integrations along the 
mean fiux direction, resulting from the emphasis 
on strong fields from the density-weighted line of 
sight integration. For integrations perpendicular 
to the mean fiux direction, the Zeeman measure- 
ments tend to underestimate the field strength. 
Thus, Zeeman measurements cannot generally be 
regarded as a lower limit on the field strength. 



A major shortcoming of our model is its low field 
strength, which renders its regime less magnetized than 
observed in molecular cores. While this is definitely a 
limitation, we see it also as an advantage since the CF- 
method is expected to work more reliably for stronger 
field strengths (HZMLN). Obviously, further limitations 
of this study stem from the limited resolution within the 
cores. 

Taking the emissivity proportional to the density 
(eq. [TJ in the Zeeman estimates does not account for 
the fact that the tracer species do not necessarily con- 
tribute proportionally to gas density. In that sense, the 
overestimate inferred for the Zeeman measurements from 
our models could be regarded as an upper limit. We did 
not address the problem of how to convert tracer column 
densities into H2 densities, which is very likely another 
major source of uncertainty in the CF-method. Similarly, 
we left out a discussion of the velocity tracers. 

We conclude that single magnetic field estimates in 
molecular cores should be regarded with caution because 
of the substantial scatter. The CF-method works reliably 
on average and should be combined with other types of 
estimates. 
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